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We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an 
adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code. We seek 
to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar 
surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split 
into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the 
poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in 
the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located 
at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, 
especially when the tube is introduced at low latitudes. 
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1 Introduction 

The Sun possesses a large variety of magnetic phenomena 
such as sunspots following an 11 -yr cycle, explosive flares 
or CME's. Such active regions on the solar surface are be- 
lieved to take their origin from strong toroidal fields created 
at the base of the convection zone (CZ). We thus need to 
understand the rising mechanisms of strong toroidal struc- 
tures through the CZ. Many models carried out since the 
80's relied on the assumption that toroidal flux is organized 
in the form of discrete flux tubes which will rise cohesively 
from the base of the CZ up to the solar surface (see Catta- 
neo, Brummell & Cline 2006 however for a less idealized 
idea of the topology of buoyant flux structures). The first 
models using "thin flux tube approximation" (Spruit 1981) 
showed that the initial strength of magnetic field was an im- 
portant parameter in the evolution of the tube. They also 
showed that an uncompensated magnetic curvature force 
could deflect the trajectory of the tube poleward, making the 
tube emerge at very high latitudes, contrary to what is ob- 
served in the Sun (Spruit & van Ballegooijen 1982; Moreno- 
Insertis, Schiissler & Ferriz Mas 1992). These results were 
also obtained in axisymmetric simulations of flux tubes in 
a rotating background (Choudhuri & Gilman 1987). More 
sophisticated multidimensional models in cartesian geom- 
etry were then developed (see review of Fan, 2004). These 
calculations showed that a sufficient twist is needed to main- 
tain the coherence of the tube during its rise (e.g. Emonet & 
Moreno-Insertis 1998). We here present the first attempt to 
study the 3-D full MHD evolution of toroidal flux tubes in 
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spherical geometry, focusing on the effects of the twist of 
the field lines, of rotation and of the initial latitude of the 
tube on its emergence. 

The paper is organized in the following manner. In sec- 
tion 2, we present the initialization of the model; in section 
3 we discuss the effects of the twist of the field lines, of a 
background rotation and of the initial latitude on the dynam- 
ical evolution of the flux tube and we conclude in section 4. 

2 The model: initial conditions and 
background state 

We use the ASH code (anelastic spherical harmonic, see 
Brun, Miesch & Toomre 2004 for details) to solve in three 
dimensions the MHD equations in a convective spherical 
shell. ASH uses the anelastic approximation that filters sou- 
nd waves without suppressing the effects of density stratifi- 
cation. Solar values are used for the rotation rate and lumi- 
nosity, and the initial stratification is obtained from a one- 
dimensional solar structure model. All the simulations are 
computed in a shell with r 6 [0.7, 0.96]i?0 with the resolu- 
tion N r = N g = Ay 2 = 256. 

To compute the model, we introduce at the starting time 
a torus of magnetic field in entropy equilibrium with the sur- 
rounding medium at the base of the computational domain 
and we let this MHD-simulation evolve. In this first paper, 
we will not address how such coherent idealized magnetic 
flux tubes are created within the Sun. This regular axisym- 
metric magnetic structure is embedded in an unmagnetized 
stratified medium. We choose for simplicity to keep the tube 
axisymmetric. We could have introduced an inhomogenous 

© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



790 



L. Jouve & A. S.Brun: Rising flux tubes 



tube with respect to longitude in our 3-D simulation, to 
study for example f2-loops emergence. However, our main 
goal here is to have a reference case to which we will com- 
pare our fully convective non isentropic simulations in a fol- 
low up paper (Jouve & Brun 2008). In order to keep a di- 
vergenceless magnetic field, we use a toroidal-poloidal de- 
composition, 



B = V x Vx(Ce r ) + Vx(Ae r 



(1) 



the expressions used for the potentials A and C for the flux 
tubes are: 



A = —Bq r exp 



r-R t 



1 + tanh 2 



(2) 



a/R t 



C 



-5 y q exp 



(3) 



1 + tanh 2 



a/R t 

where Bq is a measure of the initial field strength, a is the 
tube radius, (R t ,0 t ) is the position of the tube center and q 
is the twist parameter. 

The entropy gradient of the background state is set to 
zero so that we do not trigger the convection instability in 
this study. The density contrast between the top and the bot- 
tom of the CZ is here about 40 and the tube radius is set 
to 2 x 10 9 cm, about 0.4 times the pressure scale height at 
the base of the CZ. We compute an untwisted case (the ini- 
tial field is exclusively oriented in the direction of the tube, 
i.e. q=Q), a twisted case (with a twist above the threshold of 
Eq. 4), cases with tubes located at three different latitudes 
(namely at 15°, 45° and 60°) in a rigidly rotating shell and 
we will compare them with a non-rotating computation. 

3 Dynamical evolution of the tube 

Since the tube is introduced in entropy and pressure equi- 
libria it possesses a lower density than the surrounding me- 
dium. Consequently, the tube is buoyant and begins to rise. 

3.1 Untwisted vs twisted tubes 

Figure 1 shows the results of the evolution, in a shell rotat- 
ing at the rate = 2.6 x 10 -6 rad.s -1 , of a flux tube with 
an initial intensity of 1.8 x 10 5 G and located at a latitude of 
45° in both the untwisted case and in the twisted case (with 
a twist above the threshold, see discussion below), at two 
instants of the dynamical evolution. 

We clearly see in the case without twist (top row of Fig. 
1), the formation of two counter vortices, splitting the tube 
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Fig. 1 The four panels are snapshots of for the untwisted 
case (two upper panels) and the twisted case (two lower panels), 
cut at a specific longitude, shown in a northern quadrant of the 
sphere, after 2 days (left panels) and 7 days (right panels) of evo- 
lution. 



in two parts while it rises. As shown in Emonet & Moreno- 
Insertis (1998), this splitting can be understood by study- 
ing the equation for azimuthal vorticity. If we do not have 
any twist, the projection of the Lorentz force on the trans- 
verse plane, which constitutes a sink of vorticity, vanishes. 
Consequently, no processes can counteract the generation 
of vorticity due to the gravitational torque which results in 
the formation of two counter vortices as the tube rises. In 
the twisted case (lower row of Fig. 1), with a twist above 
a certain threshold, magnetic tension prevents vorticity to 
be created in the main body of the tube by counteracting 
the deformation of the field lines: the tube remains coherent 
during its rise. 

By setting the gravitational torque to be equal to the pro- 
jection of the Lorentz force in the equation for the azimuthal 
vorticity, we can determine the threshold above which the 
twist of the field lines can counteract the creation of two 
counter vortices inside the tube. Defining the pitch angle ip 
as the angle between the direction of the vector magnetic 
field and the longitudinal direction, we derive: 



simp 




■ 



where H p is the pressure scale height at the base of the CZ, 
Ap/p is the density deficit inside the tube compared to the 
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background stratification divided by the background density 
at the tube center and j3 is the plasma-/? associated with the 
tube. In our case, the threshold value is equal to 0.32 (cor- 
responding to a pitch angle of 18.7°). In the twisted case, 
we then use for sinV> a value of 0.52 (corresponding to a 
pitch angle of 31.3°), i.e. well above the threshold, so that 
the tube is able to rise cohesively through the entire CZ. 

This amount of twist can be related to the winding de- 
gree n of the field lines by the following formula: 

7ri? t sin 9 t 



2a 



■ tan tp 



(5) 



which gives in our case a winding degree of about 17. At 
the threshold, this would correspond to a winding degree of 
about 9. 

3.2 Effects of rotation 

To investigate the influence of rotation in our simulations, 
we compute a case where we remove the reference frame 
rotation. The initial latitude of the magnetic torus is still 
45°and the tube is still introduced at the base of the convec- 
tion zone, i.e. at about 0.75i?Q. Figure 2 shows the behavior 
of the tube in the non-rotating case (left panels) and in the 
rotating case after about 7 days of evolution. 



the rotating case, a retrograde zonal flow is created inside 
the tube which induces a Coriolis force directed towards the 
Sun's rotation axis which acts to deflect the trajectory of the 
tube poleward. Thus, we note that the deviation to the radial 
trajectory in this case is even more pronounced, reaching the 
value of about 3.8°. Another striking difference between the 
two configurations resides in the rise velocity. Indeed, look- 
ing at Fig. 2, we see that at about the same instant, the tube 
in the rotating case has not risen as high in the CZ as the 
tube in the non rotating case. 



Rise velocity 
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Fig. 2 Cut at a constant longitude of after about 7 days of 
evolution of a flux tube initially located at the base of the CZ at 
9 — 45° in the case of a non-rotating shell (left panel) and in a 
rotating case (right panel). The dashed lines represent the purely 
radial trajectory in each case. 



The most striking point is that in the non-rotating case, 
the tube remains much more symmetric with respect to the 
trajectory of the apex in a constant longitude plane. On the 
contrary, in the rotating case, we clearly see the distortion 
of the tube caused by the Coriolis force, the flux tends to be 
concentrated in the lobe which is closer to the rotation axis. 

Figure 2 also shows the poleward slip of the flux tube 
characterized by the deviation to the purely radial trajectory 
represented by the dashed line. In the non-rotating case, the 
latitudinal component of the magnetic curvature acts to drag 
the tube poleward as it cannot be compensated by any equa- 
torward force. The deviation angle to the radial direction is 
determined by comparing the latitude of the maximum of 

to its initial latitude. It is equal to 1.7° in this case. In 



Fig. 3 Rise velocity of the flux tube in the rotating and non- 
rotating cases. 



By following the position of the maximum of with 
respect to time, we are able to determine the profile of the 
rise velocity of the tube. Figure 3 shows the evolution of the 
rise velocity of the flux tube in both cases. We note that the 
accelerating phase is similar, the maximum speed is about 
the same (about 220 m.s" 1 ) but the decelerating phase is 
significantly modified in the rotating case. Indeed, the ra- 
dial component of the centrifugal force decreases the tube 
velocity so that after 6 days of evolution, the rise velocity of 
the tube in the non-rotating case is about 1.5 times that of 
the tube in the rotating case. This can be explained by look- 
ing at the equation of evolution of the radial velocity in the 
simplified model of a thin axisymmetric flux tube in spheri- 
cal geometry (Moreno-Insertis et al. 1992) which reads (as- 
suming that the rotation rate 57 of the tube is equal to the 
rotation rate of the surrounding medium and neglecting the 
advection terms): 



~dt 



<9 



Airrp 



2 sin 6^ft(6) 



The first term on the r.h.s of the equation is the buoyancy 
term which is modified by an extra term coming from the 
rotation and which has the effect of limiting the efficiency 
of buoyancy, thus resulting in a slower emergence in the 
rotating case. 
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3.3 Influence of the initial latitude 



4 Conclusion 



Since our computations are made in a spherical shell, we 
are able to check the influence of spherical geometry on the 
dynamical evolution of the tube. We can test for instance 
the influence of the initial latitude of the tube on its emer- 
gence. Figure 4 shows the topology of in a constant lon- 
gitude plane for flux tubes initially located at different lat- 
itudes. The snapshots are chosen such that the tubes have 
reached the same height in the CZ. We note that the tube at 
the latitude of 15° is delayed by almost one day in compar- 
ison to the tube at the latitude of 60°, we thus conclude that 
tubes located at higher latitudes, in a rotating background, 
are faster. In the non-rotating case, the rise velocity is the 
same for both tubes. It is again the modified buoyancy term 
of Eq.6 that is responsible for the lower rise velocity of the 
tube at low latitudes, since an increase in 6 causes an in- 
crease of r sin 2 6 ft 2 and thus a decrease of the efficiency 
of buoyancy to drive the tube upward. 




time=7.92 



Fig. 4 Same format as Fig. 2 for a flux tube initially located at 
the base of the CZ at a latitude of 15° and 60°, in a rotating shell. 



Figure 4 indicates that the initial latitude of the flux tube 
has also an influence on the deviation to the purely radial 
trajectory due to the poleward slip we mentioned before. In 
the case of a high latitude tube, the deviation reaches the 
value of 3.1° whereas in the low latitude case, the poleward 
slip is small, reaching only the value of 1.5°. This result 
is also true for tubes introduced in a non-rotating shell and 
can be understood by looking at the very simplified model 
of an axisymmetric thin flux tube rising in a non-rotating 
shell. As Moreno-Insertis et al. (1992) indicate, we can un- 
derstand the poleward drift in writing the equation for the 
^-component of the velocity in the non-rotating case, ne- 
glecting the advection terms: 



dvg 
dt 



.3 

Anrp 



■ cot( 



(7) 



This equation indicates that the acceleration in the 6- 
direction is proportional to cot 6 which is a decreasing func- 
tion of 8 between and ir/2. As 9 is here the colatitude, the 
acceleration at higher latitudes is thus more rapidly active 
than at low latitudes and as a consequence, the poleward 
drift is much more visible for a flux tube originally located 
at a latitude of 60°. 



We have carried out the first study of the 3-D full MHD 
evolution in a spherical shell of a buoyant magnetic flux 
rope in an adiabatically stratified convective zone, with spe- 
cial emphasis on the effects of rotation and of the initial 
latitude of the tube. We confirm several results obtained 
in different configurations such as in cartesian coordinates 
(e.g. Emonet & Moreno-Insertis 1998) or within the thin 
flux tube approximation (e.g. Caligari, Moreno-Insertis & 
Schiissler 1995). To be able to rise cohesively up to the solar 
surface, we confirm that a sufficient amount of twist of the 
field lines is needed. Moreover, we see in our simulations 
that twisted flux tubes can be strongly influenced by a back- 
ground rotation which has two major effects: decreasing the 
rise velocity through the spherical shell and distorting the 
flux tube such that more flux gets concentrated closer to the 
rotation axis. In all these simulations, a poleward drift of 
the tube is observed due to the non-compensated curvature 
force linked to the magnetic tension existing in the tube. 
Nevertheless, the spherical nature of these calculations en- 
ables to show that this poleward drift is minimized in the 
case where the tubes originate at low latitudes. This leads to 
consider that these low latitudes flux tubes could emerge at 
the surface in the observed solar activity belt. In this present 
work where flux tubes evolve in an isentropic background, 
3-D effects do not play any role as our tubes remain ax- 
isymmteric during their rise. On the contrary, interactions 
with fully developed convection will lead to a modulation 
of the magnetic field in the azimuthal direction and thus to 
a creation of particularly active longitudes at the solar sur- 
face. The model has thus now to be improved especially 
by computing a fully developed convective case which will 
enable us to study the non-linear interactions of these flux 
tubes with self-consistently generated differential rotation, 
meridional flows and convective plumes. Such efforts are 
under way (Jouve & Brun 2006, 2008). 
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